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Abstract A high fidelity model is developed for an 
elastic string pendulum, one end of which is attached 
to a rigid body while the other end is attached to an 
inertially fixed reel mechanism which allows the un- 
stretched length of the string to be dynamically varied. 
The string is assumed to have distributed mass and 
elasticity that permits axial deformations. The rigid 
body is attached to the string at an arbitrary point, and 
the resulting string pendulum system exhibits nontriv- 
ial coupling between the elastic wave propagation in the 
string and the rigid body dynamics. Variational meth- 
ods are used to develop coupled ordinary and partial 
differential equations of motion. Computational meth- 
ods, referred to as Lie group variational integrators, are 
then developed, based on a finite element approxima- 
tion and the use of variational methods in a discrete- 
time setting to obtain discrete-time equations of mo- 
tion. This approach preserves the geometry of the con- 
figurations, and leads to accurate and efficient algo- 
rithms that have guaranteed accuracy properties that 
make them suitable for many dynamic simulations, es- 
pecially over long simulation times. Numerical results 
are presented for typical examples involving a constant 
length string, string deployment, and string retrieval. 
These demonstrate the complicated dynamics that arise 
in a string pendulum from the interaction of the rigid 
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body motion, elastic wave dynamics in the string, and 
the disturbances introduced by the reeling mechanism. 
Such interactions are dynamically important in many 
engineering problems, but tend be obscured in lower 
fidelity models. 

Keywords Lagrangian mechanics • geometric inte- 
grator • variational integrator • string pendulum • reel 
mechanism • rigid body 



1 Introduction 

The dynamics of a body connected to a string appear 
in several engineering problems such as cable cranes, 
towed underwater vehicles, and tethered spacecraft. Sev- 
eral types of analytical and numerical models have been 
developed. Lumped mass models, where the string is 
spatially discretized into connected point masses, are 
developed in [UEIE]. Finite difference methods in both 
the spatial domain and the time domain are applied 
in ^l5j . Finite element discretizations of the weak form 
of the equations of motion are used in [5l[6]. Variable- 
length string models also have been developed: a vari- 
able length string is modeled based on a continuous 
plastic impact assumption in [7l|8], and a reel mecha- 
nism is considered in [QlfTOj. But, the reel mechanisms 
developed in those papers are problematic. In [9], the 
deployed portion of the string is assumed to move along 
a fixed line. The dynamic model of reeling developed 
in [lOj is erroneous (this will be discussed further in 
Section [2]). Instead of a point mass, a rigid body model 
is considered in [8j , but this paper does not provide any 
computational results. Analytical and numerical mod- 
els of a rigid body connected to an elastic string appear 

in yjj. 
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The goal of this paper is to develop an analytical 
model and a numerical algorithm that can be used for 
simulation of an an elastic string attached to a rigid 
body and an inertially fixed reel mechanism, acting un- 
der a constant gravitational potential. The string has 
distributed mass; it can move in a three-dimensional 
space while deforming axially; the rigid body attached 
to the string can translate and rotate. We assume that 
the point where the string is attached to the rigid body 
is displaced from the center of mass of the rigid body so 
that there exist nonlinear coupling effects between the 
string deformation dynamics and the rigid body dy- 
namics. The reel mechanism is an inert ially-fixed sys- 
tem consisting of a cylindrical reel on which the string 
winds and unwinds and a guide way that acts as a pivot 
for the deployed portion of the string. The portion of 
the string on the reel mechanism is assumed to be inex- 
tensible. The combined system of the string, the rigid 
body and the reel mechanism provides a realistic and 
accurate dynamic model of cable cranes and towing sys- 
tems. 

In this paper, we first show that the governing equa- 
tions of motion of the presented string pendulum can 
be developed according to Hamilton's variational prin- 
ciple. The configuration manifold of the string pendu- 
lum is expressed as the product of the real space R 
representing the configuration of the reel mechanism, 
the space of connected curve segments on describing 
the deployed portion of the string, and the special or- 
thogonal group SO (3) defining the attitude of the rigid 
body [12] . The variational principle is carefully applied 
to respect the geometry of the Lie group configuration 
manifold. We incorporate an additional modification 
term, referred to as the Carnot energy loss term p!3j . 
in the variational principle to take account of the fact 
that the portion of the string in the reel mechanism 
is inextensible. The resulting Euler-Lagrange equations 
are expressed as coupled partial and ordinary differen- 
tial equations. 

The second part of this paper deals with a geomet- 
ric numerical integrator for the model we presented for 
the string pendulum. Geometric numerical integration 
is concerned with developing numerical integrators that 
preserve geometric features of a system, such as in- 
variants, symmetry, and reversibility [14]. The string 
pendulum is a Lagrangian/Hamiltonian system evolv- 
ing on a Lie group. When numerically simulating such 
systems, it is critical to preserve both the symplectic 
property of Hamiltonian flows and the Lie group struc- 
ture for numerical accuracy and efficiency [15]. A geo- 
metric numerical integrator, referred to as a Lie group 
variational integrator, has been developed for a Hamil- 
tonian system on an arbitrary Lie group and it has been 



applied to several multibody systems ranging from bi- 
nary asteroids to articulated rigid bodies and magnetic 
systems in [16, 17j. 

This paper develops a Lie group variational inte- 
grator for the proposed string pendulum model. This 
extends the results presented in [17J by incorporating 
deformation of the string using a finite element model 
and by including a discrete-time Carnot energy loss 
term. The proposed geometric numerical integrator pre- 
serves the symplectic structure and momentum maps, 
and exhibits desirable energy conservation properties. 
It also respects the Lie group structure of the configu- 
ration manifold, and avoids the singularities and com- 
putational complexities associated with the use of lo- 
cal coordinates, explicit constraints or projection. As a 
result, this computational approach can represent arbi- 
trary translations and rotations of the rigid body and 
large deformations of the string. 

In summary, this paper develops an analytical model 
and a geometric numerical integrator for a string pen- 
dulum attached to a rigid body and a reel mechanism. 
These provide a realistic mathematical model for teth- 
ered systems and a reliable numerical simulation tool 
that characterizes the nonlinear coupling between the 
string dynamics, the rigid body dynamics, and the reel 
mechanism accurately. The proposed high-fidelity com- 
putational framework can be naturally extended to for- 
mulating and solving control problems associated with 
string deployment, retrieval, and vibration suppression 
as in [Is] . 

This paper is organized as follows. A string pendu- 
lum is described and the corresponding Euler-Lagrange 
equations are presented in Section [21 A Lie group vari- 
ational integrator is derived in Section [3l followed by 
numerical examples and conclusions in Section [4] and [51 



2 Euler-Lagrange Equations 

2.1 String Pendulum Model 

Consider a string that is composed of mass elements 
distributed along a curve. The string mass elements 
can translate in a three-dimensional space, and it is 
deformable along its axial direction. The bending stiff- 
ness of the string is not considered as the diameter of 
the string is assumed to be negligible compared to its 
length. The free end of the string is attached to a rigid 
body that can translate and rotate, and the point where 
the string is attached to the rigid body is displaced from 
the center of mass of the rigid body so that the dynam- 
ics of the rigid body is coupled to the string deforma- 
tions and displacements. The other end of the string 
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is connected to an inertially-fixed reel mechanism com- 
posed of a drum and a guide way. The string is wound 
around the drum at a constant radius, and the string 
on the drum and in the guide way is assumed to be in- 
extensible. A control moment is applied at the rotating 
drum. This system of the string, the rigid body, and the 
reel mechanism, acting under a constant gravitational 
potential, is referred to as a string pendulum. This is 
illustrated in Figure [TJ 

We choose an inertially fixed reference frame and a 
body- fixed frame. The origin of the body- fixed frame 
is located at the end of the string where the string is 
attached to the rigid body, and it is fixed to the rigid 
body. Since the string is extensible, we need to distin- 
guish between the arc length for the stretched deformed 
configuration and the arc length for the unstretched ref- 
erence configuration. Define 

Vd G the location of the origin of the axis 

of the drum 
d G M the radius of the drum 

Id = f^dd"^ ^ ^ the rotational inertia of the drum 

for 

6 G M the length from the drum to the 

guide way 

G M the control moment applied at drum 

L G M the total unstretched length of the 

string 

/i G M the mass of the string per the unit 

unstretched length 

O the point at which the string is at- 

tached to the drum 

s G [0, L] the unstretched arc length of the 

string between the point O and a 
material point P on the string 
t) G M the stretched arc length to a mate- 

rial point P 

Sp{t) G [6, L] the arc length of the string between 
the point O and the material point 
on the string located at the guide 
way entrance 

r(s, t) G the deformed location of a material 

point P 

0{s) gR 0= {{sp-b)-s)/dfoTs e [0,Sp-b] 

M G M the mass of the rigid body 

J G the inertia matrix of the rigid body 

with respect to the body fixed frame 
Pc G the vector from the origin of the 

body fixed frame to the center of 
mass of the rigid body represented 
in the body fixed frame 
R G SO (3) the rotation matrix from the body 

fixed frame to the reference frame 



i? G 



the angular velocity of the rigid body 
represented in the body fixed frame 



The configuration of the string on the drum and in 
the guide way is completely determined by the variable 
Sp{t), since the string there is inextensible. The con- 
figurations of the deployed portion of the string and 
the rigid body are described by the curve r(s, t) for 
s G and the rotation matrix R G S0(3), respec- 

tively, where the special orthogonal group is SO (3) = 
{R G R^""^ I R^R = I,det[R] = 1}. Therefore, the con- 
figuration manifold of the string pendulum is the prod- 
uct of the real space R, the space of connected curves 
on R^, and the special orthogonal group SO (3). 

The attitude kinematics equation of the rigid body 
is given by 



R — Rf2^ 



(1) 



where the hat map ^ : R^ ^ 50 (3) is defined by the 
condition that xy = x x y for any x^y G R^. Since x 
is a 3 X 3 skew-symmetric matrix, we have x 



T _ 



-X. 



The inverse map of the hat map is referred to as the 
vee map: {-y : 5o(3) R^. 



2.2 Lagrangian 

We develop Euler-Lagrange equations for the string pen- 
dulum according to Hamilton's variational principle. 
The Lagrangian of the string pendulum is derived, and 
the corresponding action integral is defined. Due to the 
unique dynamic characteristics of the string pendulum, 
the variation of the action integral should be carefully 
developed: (i) since the unstretched length of the de- 
ployed portion of the string is not fixed, when deriv- 
ing the variation of the corresponding part of the ac- 
tion integral, we need to apply Green's theorem; (ii) 
since the attitude of the rigid body is represented in 
the special orthogonal group, the variation of rotation 
matrices are carefully expressed by using the exponen- 
tial map pT| fT6]: (iii) since the portion of the string 
on the guide way and the drum is inextensible, the ve- 
locity of the string is not continuous at the guide way 
entrance. To take account of the effect of this velocity 
discontinuity, an additional modification term, referred 
to as a Carnot energy loss term is incorporated p!3] . 
Then, Euler-Lagrange equations are derived according 
to Hamilton's principle, and they are expressed as cou- 
pled ordinary and partial differential equations. 

Lagrangian The total kinetic energy is composed of the 
kinetic energy of the portion of the string on the drum 
and the guide way T^, the kinetic energy of the deployed 
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global reference frame 




(a) Reference configuration 
Fig. 1 String Pendulum Model 



s = 




(b) Deformed configuration 



portion of the string Tg, and the kinetic energy of the 
rigid body T^. The kinetic energy Tr can be written as 



in the reference frame is given by r{L) -\- Rp. Then, the 
kinetic energy of the rigid body is given by 



T — 



\jir{-s)-f{-s)d-s^\lde{S)f, 



where the dot represents the partial derivative with re- 
spect to time. Here, the dependency of variables on time 
t is omitted for simplicity, i.e. r(s) = r{s^t). The veloc- 
ity of the string in the reel mechanism is equal to Sp as 
the string is inextensible. From the definitions, we have 
^(0) = Sp/d, and Id = i^dd'^- Then, the kinetic energy 
Tr can be written as 



1 



Tr = ^{pSp^Kd)sl' 



(2) 



The kinetic energy of the deployed portion of the string 
is given by 



J St) 



(3) 



Let p G be the vector from the free end of the string 
r(L) to a mass element of the rigid body, expressed in 
the body fixed frame. The location of the mass element 
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-\\r{L) + Rf2pYdm 



body 



= ^Mr(L) . r{L) + Mr{L) • RQpc + ^i? • Ji7, (4) 

where J = j —p^ dm is the inertia matrix of the rigid 
body in the body fixed frame. 

Now we obtain expressions for the potential energy 
of each part. The gravitational potential energy of the 
portion of the string on the drum and the guide way is 
given by 



Vr = — • 63 — dsinO) ds 

Jo 

= -Jig {{sp -d)rd-es^ d'^{cos{{sp - h)/d) - 1)} . 

(5) 

The strain of the string at a material point located at 
r(s) is given by 



e = lim = s (s) - 1, 

As^O As V / ' 
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where ( denote the partial derivative with respect to 
s. The tangent vector at the material point is given by 



dr{s) dr{s) ds r\s) 



ds 



ds ds{s) s^{s)' 



Since this tangent vector has unit length, we have s'{s) = 
||r'(5)||. Therefore, the strain of the string is given by 
e = ||r'(s)|| — 1. The potential energy of the deployed 
portion of the string is composed of the elastic potential 
and the gravitational potential energy: 



Vs= I -EA{\y{s)\\ -If -Tigris) -esds, (6) 



where E and A denote the Young's modulus and the 
cross sectional area of the string, respectively. The grav- 
itational potential of the rigid body is given by 



Vb = -Mg{r{L)^Rpc)-e3. 



(7) 



In summary, the Lagrangian of the string pendulum 
is given by 

L={Tr- Vr) + (T, - Vs) + (Tfc - Vk) = Lr + L, + Lb. 

(8) 



2.3 Variational Approach 

Action Integral The action integral is defined by 







/ Lr 

J to 



■Ls^Li,dt = (3r^^s^^b' 



(9) 



We find expressions for the variation of each term of 
the action integral. 



Variation of 0^ From ([2]) and ([5]), the variation of 0^ 
is given by 



- {/iSp + Kd)sp - + iJ^g (rd ■ es) 



— JIgdsm{{sp — h)/d) j> 5sp dt, 
where we used integration by parts. 



(10) 



Variation of (3s From ([3|), (|6|), ([9]), the second term of 
the action integral 05 is a double integral on G 
[to^tf] X [sp{t), L]. Since the variable Sp{t) is dependent 
on the time t, the variation of 0^ should take into ac- 
count the variation of Sp{t): 



S(5s 



/ / JIr{s) • Sr{s) 

-'to Jsr,(t) 

r\s) ■ Sr\s) -{-Jlges ■ Sr{s) dsdt 



EA 



to 'J Sp{t) 



' ^^Tlf{s+)-f{s+)-\EA{\\r'{s+)\\-lf 



+ iigr{sp) ■ 63 \ 5sp dt 



(11) 



where r(s+) represents the material point of the string 
located just outside the guide way. 

Now we focus on the first term of (pT|) . Here, we 
cannot apply integration by parts at time t, since the 
order of the integrals in (fTT]) cannot be interchanged due 
to the time dependence in the variable Sp{t). Instead, 
we use Green's theorem, 



r{s) • Sr{s) ds 



to J Sp{t) 



d_ 

'dt 



{r{s) ■ Sr{s)) dsdt, 



(12) 



where represents the counterclockwise line integral 
on the boundary B of the region [to^tf] x [sp{t), L]. The 
boundary B is composed of four lines: (t = to^s G 
[sp{to),L]), {t = tf,s e [sp{tf),L]), {t e [to,tf],s = 
Sp{t)), and {t G [to,t/],s = L). For the first two lines, 
Sr{s) = since t = to,tf. For the last line, ds = since 
s is fixed. Thus, parameterizing the third line by t, we 
obtain 



r{s) - Sr{s) ds = / f{sp{t))-5r{sp{t))sp{t)dt. 
B Jto 



Substituting this into ([12]) and rearranging, the first 
term of (pT]) is given by 



r{s) ' Sr{s) 



to 



-r{s) ' 6r{s) ds + r{sp) • 6r{sp) Sp 



dt. 



(13) 



Substituting this into (pTj) . and using integration by 
parts with respect to s for the second term of (fTT]) . the 
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variation of 0^ can be written as 

rtf pL 
I to 



6(5, 



ds dt 



= f ' f {-Ms) + F\s) + 71^63} • Sr{s) 

J to J Sp 

f^' { - \w{st) ■ r{st) + \EA{\\r'{s^)\\ - if 



F{L)-5r{L)+F{sp)-5r{s+)dt, 



where F{s) = " ||^,^Jj|| r'{s) represents the tension 
of the string. 

We simphfy this using the boundary condition at 
the guide way. The location of the guide way entrance is 
given by = r{sp{t), t). Since the location is inertially 
fixed, we have Svp = Sr{s^) + r\s^)Ssp = 0, and f-p = 
r(s+) + r\s^)sp = 0. Substituting these, we obtain 

8(5 s = / {-Wi'^) + F\s) + JIges} • Sr{s) dsdt 

J to J Sp{t) 

+ £' \]^Tl\\r\s;)f sl + \EA{\\r\s;)\\-lf 



figr{sp) ■ 63 1 Sspdt 



r F{sp) • Ssp + F{L) • Sr{L) dt. 

J to 



(14) 



Variation of 0^ From (|4]), ([7|), the variation of 0^ is 
given by 

S(5b= r {Mr{L)^MRQpc}-5f{L) 

J to 

+ { Ji? + MpcR^r{L)} ■ Sf2 + Mr{L) ■ 5RQpc 

+ Mge^ • 5r{L) + Mge^, ■ SRpc dt. (15) 

The attitude of the rigid body is represented by the 
rotation matrix R G S0(3). Therefore, the variation 
of the rotation matrix should be consistent with the 
geometry of the special orthogonal group. In p!7|[T6], it 
is expressed in terms of the exponential map as 



SR = 



R' = 



R exp ef] = Rfj 



(16) 



for ?7 G M^. The key idea is expressing the variation 
of a Lie group element in terms of a Lie algebra ele- 
ment. This is desirable since the Lie algebra 50 (3) of 
the special orthogonal group, represented by 3 x 3 skew 
symmetric matrices, is isomorphic as a Lie algebra to 
using the vee map. As a result, the variation of the 
three-dimenstional rotation matrix R is expressed in 



terms of a vector 77 G R^. We can directly show that 
(dni) satisfies S{R^R) = SR^R + R^SR = -7) + 7) = 0. 
The corresponding variation of the angular velocity is 
obtained from the kinematics equation ([1]): 



Sf2 = 



de 



{R'fR' = {r]^f2xr])''. 



(17) 



Substituting (p!6|) . (pTj) into (p!5]) . and using integra- 
tion by parts for 77, we obtain 

S(5b= r M{-r{L) - RQ'^pc - RQpc + ge^^} ■ r{L) 

Jto 

+ {-Jf2 - MpcR^r{L) - QJQ + MgpcR^es} -rjdt. 

(18) 



T8l). the variation of 



Variation of (5 From ([TOl), ([T 
the action integral is given by 

6(5 = 6(5r^S(5s^S(5b. 



(19) 



Variational Principle with Discontinuity Let = r{sp{t)^t) 
be the location of the pivot in the reference frame. Since 
it is fixed, we have rp = r{sp^t) + r'{sp^t)sp = 0. Let 
r(5~), and r(s+) be the material point of the string just 
inside the guide way, and the material point just outside 
the guide way, respectively. Since the string is inextensi- 
ble inside the guide way, ||r'(s~)|| = 1. Since the string 
is extensible outside the guide way, ||r'(s+)|| = 1 + e+, 
where e+ represents the strain of the string just outside 
the guide way. Using these, the speeds of the string at 
those points are given by 

\\r{s;)\\ = \\-r'{s;)sp\\ = \sp\ 
\\r{s;)\\ = \\-r'{.s+)sp\\ = {l + e+)\sp\. 

Therefore, the speed of the string changes instanta- 



neously by the amount 



at the guide way. 



Due to this velocity and strain discontinuity, the 
variation of the action integral is not equal to the neg- 
ative of the virtual work done by the external control 
moment u at the drum. In order to derive equations of 
motion using Hamilton's principle, an additional term 
Q, referred to as Carnot energy loss term should be in- 
troduced [ISIIH]. The resulting variational principle is 
given by 



5(5 



/ {Q 

^to 



u/d)Ssp dt = 0. 



(20) 



The corresponding time rate of change of the total en- 
ergy is given by E = {Q -\- u/d)sp^ where the first term 
Qsp represents the energy dissipation rate due to the 
velocity and strain discontinuity. 
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Consider the infinitesimal mass element dm = JlSpdt 
located just outside the guide way. Without loss of gen- 
erality, we assume that Sp > (retrieval case). The 
motion of this mass element moving with the velocity 
(1 + €~^)sp can be considered as a plastic impact into 
the the portion of the string on the guide way moving 
with velocity Sp. The corresponding energy dissipation 
rate is given by 



(21) 



(See [inilH]). Dividing both side by Sp^ we obtain the 
expression for the Carnot energy loss term Q. 



this case, the total energy and the total angular mo- 
mentum about the gravity direction 63 are conserved: 

E={Tr^ Vr) + {Ts + Vs) + (T5 + H), 

/ fir{s) X r{s) ds 

. J Sr, 



Mr{L) X (r(L) + ROpc) - Mr{L) x Rp^ + Ji7 



If we choose pc = 0^ then the rotational dynamics of 
the rigid body ([25]) is decoupled from the other equa- 
tions. In this case, ([22|) - ([24|) describe the dynamics of 
an elastic string attached to a point mass and a reel 
mechanism. 



2.4 Euler-Lagrange Equations 

Substituting (p!9|) . ([2T]) into ([2Q|) . we obtain Euler-Lagrange 
equations for the string pendulum: 

-{jtSp + i^d)sp + Jig {vd - Tp) • 63 - JIgdsm{{sp - h)/d) 

- F{s;) • +7i(||/(4)|| - 1)^^ + ^ = 0, 

(22) 

-/Ir{s) + F\s) + Jlges = 0, se[sp,L], (23) 

-Mr{L) - MRQ^pc - MR^pc + Mges - F{L) = 0, 

(24) 

-jn - MpcR^r{L) - QJQ + MgpcR^es = 0, (25) 

where F{s) = EA " r'{s) is the tension of the 

string. These are coupled ordinary and partial differen- 
tial equations. The motion of the reel mechanism and 
the deployed portion of the string are described by ([22]) 
and ([23]) , respectively. The translational and rotational 
dynamics of the rigid body are determined by ([24]) . ([25]) . 
All of these equations are coupled. 

In ([22]) . the fifth term, /le+s^, represents the effect of 
the velocity discontinuity. Note that this term vanishes 
if the deployed portion of the string is also inextensi- 
ble, i.e. ||r'(5)|| = 1 or e+ = 0. A similar expression is 
developed in [10] from momentum balance, but their 
expression is erroneous. 

Special Cases Suppose that the length of the string on 
the reel mechanism is fixed, i.e. Sp{t) = Sp{to) for any 
t > to. Then, the equations of motion ([23 ]) - ([25]) describe 
the dynamics of an elastic string pendulum model with 
a fixed unstretched length, which is studied in [TT]. In 



3 Lie Group Variational Integrator 

Geometric numerical integration deals with numerical 
integration methods that preserve geometric properties 
of a dynamic system, such as invariants, symmetries, re- 
versibility, or structure of the configuration manifold [HI 
19]. The geometric structure of a dynamic system de- 
termines its qualitative dynamical behavior, and there- 
fore, the geometric structure-preserving properties of a 
geometric numerical integrator play an important role 
in the qualitatively accurate computation of long-term 
dynamics. The continuous-time Euler-Lagrange equa- 
tions developed in the previous section provide an ana- 
lytical model for a string pendulum. However, the pop- 
ular finite difference approximations or finite element 
approximations of those equations using a general pur- 
pose numerical integrator may not accurately preserve 
the geometric properties of the system [14 . 

Variational integrators provide a systematic method 
of developing geometric numerical integrators for La- 
grangian/Hamiltonian systems [20]. Discrete-time Euler- 
Lagrange equations, referred to as variational integra- 
tors, are constructed by discretizing Hamilton's princi- 
ple, rather than discretizing the continuous- time Euler- 
Lagrange equations using finite difference approxima- 
tions. This is in contrast to the conventional viewpoint 
that a numerical integrator of a dynamic system is a 
discrete approximation of its continuous-time equations 
of motion. As it is derived from a discrete analogue of 
Hamilton's principle, it preserves symplecticity and the 
momentum map, and it exhibits good total energy be- 
havior for an extremely long time period. 

On the other hand. Lie group methods conserve the 
structure of a Lie group configuration manifold as it up- 
dates a group element using the group operation [21]. 
As opposed to computational methods based on local 
coordinates, projections, or constraints, this approach 
preserves the group structure naturally without any 
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singularities associated with local coordinates or the 
additional computational overhead introduced by con- 
straints. 

These two methods have been unified to obtain a Lie 
group variational integrator for Lagrangian/Hamiltonian 
systems evolving on a Lie group [17]. This geometric in- 
tegrator preserves symplecticity and group structure of 
those systems concurrently. It has been shown that this 
property is critical for accurate and efficient simulations 
of rigid body dynamics [15j. This is particularly useful 
for dynamic simulation of a string pendulum that un- 
dergoes large displacements, deformation, and rotations 
over an exponentially long time period. 

In this section, we develop a Lie group variational 
integrator for a string pendulum. We first construct a 
discretized string pendulum model, and derive an ex- 
pression for a discrete Lagrangian, which is substituted 
into discrete-time Euler-Lagrange equations on a Lie 
group. 



3.1 Discretized String Pendulum Model 

Let > be a fixed time step. The value of variables 
dit t = to -\- kh is denoted by a subscript k. We discretize 
the deployed portion of the string using N identical line 
elements. Since the unstretched length of the deployed 
portion of the string is L — Sp^^ the unstretched length of 
each element \s Ik = — W^- Let the subscript a denote 
the variables related to the a-th element. The natural 
coordinate of the a-th element is defined by 



(26) 



for s G [spf^ + (a — l)//e, +a//e]. This varies between 
and 1 for the a-th element. Let ^o. Si be shape functions 
given by So{C) = 1 — and Si{C) = These shape 
functions are also referred to as tent functions. Define 
Qk^a to be the relative location of a string element with 
respect to the guide way entrance, i.e qk^a = ^/e,a — ^p- 
Using this finite element model, the position vector 
r(5, t) of a material point in the a-th element is approx- 
imated as follows: 



rk{s) = So{Ck,a)qk,a + Si{Ck,a)qk,a 



+1 



(27) 



Therefore, a configuration of the presented discretized 
string pendulum dit t = kh -\- to is described by gk = 
{spf^;qk,i, . . . , g/c,Ar+i; Rk), and the corresponding con- 
figuration manifold is G = R x (R3)^+i x S0(3). This is 
a Lie group where the group acts on itself by the diago- 
nal action [12] : the group action on Sp^ and Qi^k ' ' ' QN-\-i,k 
is addition, and the group action on Rk is matrix mul- 
tiplication. 



We define a discrete-time kinematics equation using 
the group action as follows. Define fk = (^^^^ ; Aqk^i^ . . ., 
Z\g/c,Ar+i; F/c) G G such that gk^i = gkfk- 

(Spfc+i; g'/c+l,!, • • • , ^/c+l,Ar+l, Rk^l) = 

i^Pk + ^Spj^;qk^i + ^^/e,i, . . . , qk,N+i + ^^/c,Ar+i; Rk^k)- 

(28) 

Therefore, fk ^ G represents the relative update be- 
tween two integration steps. This ensures that the struc- 
ture of the Lie group configuration manifold is numeri- 
cally preserved since gk is updated by fk using the right 
Lie group action of G on itself. 



3.2 Discrete Lagrangian 

A discrete Lagrangian Ld{gk^fk) : G x G ^ M is an 
approximation of the Jacobi solution of the Hamilton- 
Jacobi equation, which is given by the integral of the 
Lagrangian along the exact solution of the Euler-Lagrange 
equations over a single time step: 



nh 

iigkjk)- / L{g{t),g-\t)g{t))dt, 
Jo 



where g{t) : [0, h] ^ G satisfies Euler-Lagrange equa- 
tions with boundary conditions ^(0) = gk^ g{h) = gkfk- 
The resulting discrete-time Lagrangian system, referred 
to as a variational integrator, approximates the Euler- 
Lagrange equations to the same order of accuracy as 
the discrete Lagrangian approximates the Jacobi solu- 
tion dilEo]. 

We construct a discrete Lagrangian for the string 
pendulum using the trapezoidal rule. We first find the 
contributions of each component to the kinetic energy. 
From the given discretized string pendulum model and 
([2|) , the kinetic energy of the reel mechanism is approx- 
imated by 



(29) 



Using the chain rule, the partial derivative of rk{s) 
given by ([27|) with respect to t is given by 



rk{s) = ^\ So{Ck,a)^qk,a + Si{Ck,a)^qk,a+l 



{L - s) {qk^a - qk,a+l) 



(L-Sp^) 



As, 



] 



Substituting this into (|3|), the contribution of the a-th 
string element to the kinetic energy of the string is given 
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by 



3.3 Lie Group Variational Integrator 



Tk,a = ^MlAqk,a ■ Mk,a + ^M|Z\gfe,a+l ' Mk,a+1 
+ ^M^%^Sp, ■ Aqk,a+1 + ^Ml^ASp, ■ Aqk,a, 

(30) 

where inertia matrices are defined in Appendix lA.li 
From the attitude kinetics equations ([T]), the angular 
velocity is approximated by 

4 ^ ^Rl{Rk+i - Rk) = ^{Fk - /). 

Define a non-standard inertia matrix Jd = ^tr[J]/ — J. 
Using the trace operation and the non-standard inertia 
matrix, the last term of the kinetic energy of the rigid 
body given by ([2]) can be written in terms of as ^ i? • 
Jf2 = ^tT[QJ(if^^]' Then, the kinetic energy of the rigid 
body is given by 



^MAqu^N+i ' Rk{Fk - I)pc, 



(31) 



where we use properties of the trace operator: tr[A5] = 
iT[BA] = tT[A^B^] for any matrices A,B eR^""^. 

From ([29]), (|30]), (|3lD, the total kinetic energy of the 
discretized string pendulum is given by 



AT 



Tk = Tk^r + ^ Tk^a + Tk^b' 



(32) 



a=l 



Similarly, the total potential energy for the given dis- 
cretized string pendulum can be written as 

Vk = -fig{{sp^ - • 63 + (f{cos{{sp - h)/d) - 1)} 
^ 1 

+ ^ --l^ghes • (% + Qk^a + qk,a+l) 
a=l ^ 

lEA 



+ ^^(lkfc,a+l - qk,a\\ - hf 
I Ik 

- Mges • {qk,N+i + + Rkpc)- 



For a given discrete Lagrangian, the discrete action sum 
is given by (3d = ^k^dk- As the discrete Lagrangian 
approximates the action integral over a single discrete 
time step, the action sum approximates the action in- 
tegral. According to the discrete Lagrange-d'Alembert 
principle, the sum of the variation of the action sum 
and the discrete virtual work done by external control 
moments and constraints is equal to zero. This yields 
discrete-time forced Euler-Lagrange equations referred 
to as variational integrators ^20j. This procedure is fol- 
lowed for an arbitrary discrete Lagrangian defined on a 
Lie group configuration manifold in [17] to obtain a Lie 
group variational integrator: 



T!L 



fk-i'^fk-iFdk-i - Ad* 1 • (T*L/^ • D/^LrfJ 

k 

+ T*L^, • T>g^Ld^ + Ud^ + Qdk = 0, 



^fe+i = gkfk, 



(35) 
(36) 



where T*L : G x T*G T*G is the co-tangent lift of 
the left translation action, D / represents the derivative 
with respect to /, and Ad* : G x ^ is the co- 
Adjoint operator [T2] . 

The contribution of the external control moment 
and the Carnot energy loss term are denoted by Ud^^ 
and Qdk' They are defined to approximate the addi- 
tional term in the variational principle ([2Q|) that arises 
due to a discontinuity: 



(/e+l)/i 



(Q + u/ d)Ssp dt ^ {Qdk + Udk)Ss 



kh 

From ([2T]) . these are chosen as 

Qd, = -^mslje + EA){\\qkA - Ik?, 
= huk/d. 



(37) 
(38) 



(33) 



We substitute the expressions for the discrete La- 
grangian the Carnot energy loss term (|37|). and the 
control moment ([38]) into (|35|) and (|36|) to obtain a Lie 
group variational integrator for the discretized string 
pendulum model. This involves deriving the derivatives 
of the discrete Lagrangian and their co-tangent lift. The 
detailed procedure and the resulting expressions for Lie 
group variational integrators are summarized in the Ap- 
pendix. 



This yields the discrete-Lagrangian of the discretized 
string pendulum 



ikidk.fk) = hTk{gk,fk) - -^Vk{gk,fk) - -^Vk+i{gk, fk)' 



(34) 



Computational Properties The proposed Lie group vari- 
ational integrators have desirable computational prop- 
erties. The Lie group configuration manifold is often 
parameterized. But, the local parameterizations of the 
special orthogonal group, such as Euler angles or Ro- 
drigues parameters, have singularities. In a numerical 
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simulation of large angle maneuvers of a string pen- 
dulum, these local parameters should be successively 
switched from one type to another in order to avoid 
their singularities. They also lead to excessive complex- 
ity. Non-parametric representations such as quaternions 
also have associated difficulties: there is an ambiguity 
in representing an attitude since the group SU(2) of 
quaternions double cover S0(3). Furthermore, as the 
unit length of quaternions is not preserved in numeri- 
cal simulations, attitudes cannot be determined accu- 
rately. Sometimes, numerical solutions updated by any 
one-step integration method are projected onto the Lie 
group at each time step [22]. Such projections may de- 
stroy the desirable long-time behavior of one-step meth- 
ods, since the projection typically corrupts the numeri- 
cal results. Lie group variational integrators update the 
group elements by using a group operation. Therefore, 
the Lie group structure is naturally preserved at the 
level of machine precision, and they avoid any singular- 
ity and complexity associated with other approaches. 

Since Lie group variational integrators are constructed 
according to Hamilton's principle, their numerical tra- 
jectories preserve a symplectic form and a momentum 
map associated with any symmetry. These ensure long- 
term structural stability and avoid artificial numerical 
dissipation. These properties are difficult to achieve in 
conventional approaches based on finite difference ap- 
proximation of continuous equations of motion. 

In summary, the proposed Lie group variational in- 
tegrators for a string pendulum will be particularly use- 
ful when studying nontrivial maneuvers that combine 
large elastic deformations and large rigid motions accu- 
rately over a long time period. 

4 Numerical Examples 

We now numerically demonstrate the computational 
properties of the Lie group variational integrators de- 
veloped in the previous section. The properties of the 
reel mechanism are as follows: b = d = 0.5 m, a^^^ = 1 kg. 
The material properties of the string are chosen to rep- 
resent a rubber string [5]: /I = 0.025 kg/m, EA = 40 N, 
L = 100 m. The rigid body is chosen as an elliptic cylin- 
der with a semimajor axis 0.5 m, a semiminor axis 0.4 m, 
and a height 0.8 m. The mass and the location of the 
center of mass of the rigid body are M = 0.1kg, and 
pc = [0.3,0.2,0.4] m, respectively. 

We consider three cases: (1) dynamics of a fixed 
length string pendulum released from a horizontal con- 
figuration, (2) deployment dynamics due to gravity from 
a horizontal configuration, and (3) retrieval dynamics 
using a constant control moment. Initial conditions are 
as follows: 



Spo (m) 


qo,a (m) 


Uk (Nm) 


(1) 


90 


loei 




(2) 


99 


lo{a - l)ei 





(3) 


90 


lo{a - 1) (sin 15°ei + cos 15°e3) 


2.09 



For all cases, we choose Sp — Om/s, (70, a — 

Om/s 

for a G {1,A^}, qo,N+i = 0.5e2 m/s, Rq = /, i7o = 
Or ad/sec. The deployed portion of the string is dis- 
cretized by = 20 elements, and the time step is 
h = 0.0005 second. Simulation time is T = 10, T = 8, 
and T = 10 seconds for each case, respectively. 

Energy transfer The following figures show the simu- 
lation results. The maneuver of the string pendulum is 
illustrated by snapshots, where the relative elastic po- 
tential distribution at each instant is denoted by color 
shading (the corresponding animations are available at 
http://my.fit.edu /~taeyoung). As the point where the 
string is attached to the rigid body is displaced from the 
center of mass of the rigid body, the rigid body dynam- 
ics are directly coupled to the elastic string dynamics. 
The illustrated maneuvers clearly show the nontrivial 
coupling between the strain deformation, the rigid body 
dynamics, and the reel mechanism. 

The energy exchange plots also show that there is 
significant energy transfer between the kinetic energy, 
the gravitational potential energy, and the elastic po- 



tential energy. In Figure 2(b), there is an energy ex- 
change between the kinetic energy and the gravitational 
potential energy. But, when the string is mostly stretched 
at t = 1.8 and t = 5.5 seconds, part of the kinetic en- 
ergy is transferred to the elastic potential energy and 
the rotational kinetic energy. The rigid body starts to 
tumble at t = 7 seconds. The elastic potential energy 



transfer along the string is observed in Figure [2(a) For 
the deployment case shown in Figure |3(b)l the gravi- 
tational potential energy is generally transferred to the 
kinetic energy. As the length of the deployed portion 
of the string increases, the elastic potential energy in- 
creases and the rotational kinetic energy decreases. For 
the third retrieval case, both the gravitational potential 
energy and the total energy increase due to the constant 
control moment. In addition, there is a smaller-scale 
periodic energy exchange between the elastic potential 
and the gravitational potential energy with an approx- 
imate period of 1.2 seconds. The rigid body starts tum- 
bling at t = 3 seconds. 

Conservation Properties The proposed Lie group vari- 
ational integrators exhibit excellent conservation prop- 
erties for these complicated maneuvers of the string 
pendulum. For the fixed length string dynamics, the 
total energy and the total angular momentum about 
the gravity direction should be preserved. The devia- 
tions of those quantities are shown in Figure |2(d)| where 
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the maximum deviation of the total energy is less than 
0.01% of the maximum kinetic energy, and the devia- 
tion of the angular momentum is less than 3 x 10~^% 
of its initial value. For the second deployment case, the 
total energy dissipates only due to the velocity discon- 
tinuity. Figure |3(d)| shows the difference between the 
computed total energy change and the energy dissipa- 
tion computed by the Carnot energy loss term ([2T]) . 
The difference is less than 0.0003% of the maximum 
kinetic energy, which illustrates that there is no artifi- 
cial numerical dissipation caused by the proposed Lie 
group variational integrator. The orthogonal structure 
of rotation matrices is preserved to machine precision. 
Figure [3(d)] and [4(d)] show that the orthogonality error, 
measured by || / — is less than 10~^^. 

5 Conclusions 

We have developed continuous-time equations of mo- 
tion and geometric numerical integrators, referred to as 
Lie group variational integrators, for a 3D elastic string 
pendulum attached to a rigid body and a reel mecha- 
nism. They are carefully derived while taking account 
of the length change of the deployed portion of the 
string, the Lie group configuration manifold of the rigid 
body, and the velocity discontinuity at the guide way 
entrance. The continuous-time equations of motion pro- 
vide an analytical model that is defined globally on the 
Lie group configuration manifold. The Lie group varia- 
tional integrator preserves the geometric features of the 
system, thereby yielding a reliable numerical method 
to compute the nonlinear coupling between the large 
string deformation and the nontrivial rigid body dy- 
namics accurately over a long time period. In short, 
this paper provides high fidelity analytical and compu- 
tational models for a string pendulum. 

The numerical experiments suggest that accurately 
modeling the reeling mechanism is of critical impor- 
tance in order to capture the correct dynamics, due to 
the disturbance that is introduced in the string at the 
point of contact with the reeling mechanism when the 
string is deployed or retracted. One can observe that 
this disturbance propagates down the string at a veloc- 
ity that is determined by the elastic properties of the 
string. Since the point of contact between the string 
and the rigid body does not go through the center of 
mass of the rigid body, the elastic disturbance excites 
a rotational response in the rigid body. As such, ac- 
curately modeling the reel mechanism, elastic string 
dynamics, rigid body motion, and their interactions, 
is critical for obtaining realistic predictions about how 
towed underwater vehicles and tethered spacecraft be- 
have when performing aggressive maneuvers. 



The proposed string pendulum model and compu- 
tational approach can be extended in several ways. For 
example, different types of string models can be consid- 
ered, such as an inextensible string, nonlinear elasticity, 
and bending stiffness. The reel mechanism can be gen- 
eralized by assuming that the portion of the string on 
the drum is also extensible. These results can be ex- 
tended to model tethered spacecraft in orbit, and they 
can be used to study associated optimal control prob- 
lems by adopting the discrete mechanics and optimal 
control approach [23] . 
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t 

(a) Snapshots at each 0.2 second t G [0,5] (b) Energy exchange (£^:sohd, T:sohd, Trot:dashed, V^ravity :dash- 

dotted, Vela Stic -dotted) 




t t 
(c) Stretched length of the deployed portion of the string, and the (d) Deviation of conserved quantities: total energy and the total 

second component of the angular velocity i? angular momentum about the gravity direction 



Fig. 2 Fixed length string pendulum 
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(c) Length of the deployed portion of the string (stretched:sohd, (d) Deviation of conserved quantities: the difference between the 

unstreched:dashed), and the second component of the angular computed total energy change and the energy dissipation due to 

velocity i? the velocity discontinuity, the orthogonality error of the rotation 

matrix 



Fig. 3 Deployment due to gravity 
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(a) Snapshots at each 0.5 second t G [0, 10] 





(b) Energy exchange (£^:sohd, T:sohd, 
^^ra^;ity:dash-dotted, VgZastic :dotted) 




(c) Length of the deployed portion of the string (stretched:sohd, 
unstreched: dashed), and the second component of the angular 
velocity Q 



1 2 3 4 5 6 7 

t 

(d) Orthogonality error of the rotation matrix 



Fig. 4 Retrieval using a constant control moment 
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A Development of the Lie Group Variational 
Integrator for a String Pendulum 

A.l Inertia Matrices for the Discrete Lagrangian 

The inertia matrices for the discrete Lagrangian are defined as 
fohows. 



Ml = -Tilk. Ml = Ml 
3 _ 1-7 (3^^ + 3A^ + 1 - QNa - 3a + ?>o?) 



Mia = ^7^h 

M^^-^Ul /Vf23 - l-(l + 3^-3a) 

Mk - -f^^k, ^/c,a - 



Li 

6 



{Qk,a — Qk,a-\-l)^ 



,,31 l_(2 + 3Ar-3a) , 

Mk,a = [qk,a - qk,a + l)- 

6 N 



A. 2 Derivatives of the discrete Lagrangian 

The Lie group variational integrator given by (I35|) is expressed in 
terms of the derivatives of the discrete Lagrangian and their co- 
tangent hft. Here, we describe how to compute the co-tangent hft 
and the co- Adjoint operator on the configuration manifold G = 
M X (M3)'^+^ X SO (3) without introducing the formal definition 
of those operators. 

The co-tangent lift of the left translation on a real space is 
the identity map on that real space. Using the product struc- 
ture of the configuration manifold G = M x (M^jiv+i x S0(3), 
the derivative of the discrete Lagrangian with respect to fk = 
{Aspf^;Aqk^i, . . ■ , Aqk^^^i; Fj.) G G is given by 



(39) 



Deriving the derivatives of the discrete Lagrangian with re- 
spect to Aspj^ or Aqk^a is straightforward. For example, from 
(j30|) , (|32|) , (|34|) , the derivative of the discrete Lagrangian with 
respect to Z\g/c,a for any a G {2, . . . , A^} is given by 



Similar to (jl6|) , the variation of the rotation matrix Fj^ can be 
written as SF^ = F]^C, for a vector G M^. From the definition of 
the co-tangent lift of the left translation, we have 

{~^}^Fk • ^F^LdJ • C = j^^^[-FkCkJd] + ^Rl^qk,N+i • FkCkPc 

+ ^Mge3 • RkFkCkPc- 

By repeatedly applying the following property of the trace op- 
erator, tr[AB] = tr[BA] = tr[A^B^] for any A, B e R^^^, 
the first term can be written as tT[—FkCkJd] = ^i^[—CkJdFk] = 
tr[Cfci^J Jd] = -^tT[Ck{JdFk-F^ Jd)]. Using the property of the 
hat map, x^y = —^tr[xy] for any G M^, this can be fur- 
ther written as {{Jd^k — ^k '^d)^) ' Cfc- As y • xz = zy • x for any 
x,y,z G M^, the second term can be written as F^ Aqk^N+i • 
CkPc = PcF^ R^Aqk^N+i • Ck- Using these, we obtain 

T|Li.^ • Di^^Ld^ = ^{JdFk - F^JdY + ^PcF^RlAqk,N+i 



-MgpcF^Rle3. 



(43) 



Expression for the derivatives of the discrete Lagrangian with 
respect to Rk is similarly developed. 

In summary, in addition to (|40|) and (|43|) . derivatives of the 
discrete Lagrangian are summarized as follows. 

1 1 ^ 

D^3,, Ld, = -MlAsp, + - Y.^Ml\ + Mf_,) . Aqk^a 

a— 2 

+ \Ml%Aqk,N+i - ^-Dsp^^^Vk+i, 
Mk = PSpf^ +i^d + ^P(L - SpJ, 
^Aqk,N + iLd,, = \iMl + M)Aqk,N+l + ^Ml^Aqk,N 

+ ^Mi%Asp, + ^MRk(Fk - I)pc 
h 

~ —'^qk+i,N+i^k+i^ 

N 

^-Vk = - Y.^^^k,a ■ ^qk,a 



6Nh 



^Ml^Aqk,a-i + ^MlAqk,a + -^M^ Z\gfc,,+i 



Vk. 



K,a ' K,a — ±/ Jt^K c2 

where the derivative of the potential energy is given by 



-pglke3^VVl^_^-VVl^ 



EA ||g/c,a+i - qk,a\\ - h 



\\qk,a-\-i — qk,i 



(qk,a-\-l — qk,a)- 



(40) 



(41) 
(42) 



Expressions for the other derivatives of the discrete Lagrangian 
with respect to Spj^, Aspj^, qk,a cire similarly developed and they 
are summarized later. 

Now we find the derivative of the discrete Lagrangian with 
respect to Fk- From ((3T]) . ((33]) . ((34]) . we have 

Dj.^Ld^ • 6Fk = ^tT[-6FkJd] + ^RlAqk^N+i ■ SFkPc 
+ ^Mges • Rk^FkPc. 



+ ^qk,a+i • ^qk,a+i + ^qk,a • ^qk,a+l) 

^^k,aLd^ = -^{l^3N-3a)Asp,Aqk,a+i 

- -^{5^3N -3a)Asp^Aqk,a-i 
h h 



D 



^ ^Spf^Aqk,N+l- ^^^Spf^Aqk^N 



-Dgfc,7V + lVfc - -'Dqk^N + l^k-^-l, 



T*LHfc • r>R,^Ld^ = ^-^{{Fk - I)pc)^RiAqk,N+i + 'jMgpcRies 



M 
1 
h 



MgFkPcF^Rle3, 



D 



Vk = -pgvd • es + pgdsm((spf^ - h)/d) 
1 ^ 

+ ^ 51 ^^^3 • (2rp + qk,a + qk,a+l) 
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+ -^(||^fc,a + l -gfc,af -ll), 

^lk,N+iVk = -(\lJ.lk + M)gez + Wlf,. 

The CO- Adjoint map on a real space is the identity map 
on that real space. The co- Adjoint map on SO (3) is given by 
Ad* = FkP = (FkpF^y for any p G (M^)* - 5o(3)*. Using 

the product structure of the configuration manifold, we have 
Ad* i(T*L/^-D/^Ld^)= [DAsp^Ld^;r>Aq^^^Ld^,--- , 

J k 

D4<j»+i,. Ld, ; Ad* ^ (TJ , • , Ld, )] , (44) 

k 

where 

Ad;.^^(T*Li., . r>F,LdJ = \{FkJd - JdFkY 

+ ^FkPcF^RlAqk,N+i + ^MgFkPcF^Rle3. (45) 

Discrete-time Euler- Lagrange Equations Substituting the 

derivatives of the discrete Lagrangian given in (l40T) . (1431) and the 
appendix, the co- Adjoint map given by (|44|) . and the contribu- 
tions of the external control moment and the Carnot energy loss 
term ()37|) . (|38|) into the Lie group variational integrator on an ar- 
bitrary Lie group given by p5|) . (|36|) . we obtain the discrete-time 
Euler-Lagrange equations of the string pendulum at (|46|) - (|51|) . 

Equation (|48|) is satisfied for a G {2, ... , A^}, and (|51|) is satis- 
fied for a G {1,...,A^}. For given gk = {spf^] qk,i, ■ ■ ■ qk,N+i', Rk), 
we solve (|46|) - (|50|) for the relative update = i^^Pk 5 ^Qk,i, • • -5 
Z\gfc^iV+i; ^fc)- Then, the configuration at the next step Qk^i = 
• • • , Q'fc+i,Ar+i; ^fc+i) can be obtained by (|5T]). This 
yields a discrete-time Lagrangian flow map (g^, fk) (fi'fc+i? /fc+i)? 
and this is repeated. 

Special Cases If we set Asp,^ = for any k, then the discrete- 
time Euler-Lagrange equations and Hamilton's equations provide 
a geometric numerical integrator for a string pendulum model 
with a flxed unstretched string length, studied in [TT]. If we chose 
Pc = 0, then these equations describe the dynamics of an elastic 
string attached to a point mass and a flxed pivot. 

Computational Approach These Lie group variational inte- 
grators for a string pendulum are implicit: at each time step, 
we need to solve nonlinear implicit equations to flnd the relative 
update /fc G G. Therefore, it is important to develop an effi- 
cient computational approach for these implicit equations. This 
computational method should preserve the group structure of 
fk, in particular, the orthogonal structure of the rotation matrix 
Ffc G SO (3). The key idea of the computational approach pro- 
posed in this paper is to express the rotation matrix Fj^ in terms 
of a vector G M.^ using the Cayley transformation |14) : 

Fk = {I + Ck){I - Ck)-\ (52) 

Since the rotation matrix Ff^ represents the relative attitude up- 
date between two adjacent integration steps, it converges to the 
identity matrix as the integration step h approaches zero. There- 
fore, this expression is valid for numerical simulations even though 
the Cayley transformation is a local diffeomorphism between 
and S0(3). 

Our computational approach is as follows. The implicit equa- 
tions for Fk given by (|50|) are rewritten in terms of a vector 
using (|52|) . and a relative update expressed by a vector = 



[Asp^;Aqk,i,...Aqk,N+i;ck] G R X (R3)iv+i X R3 is solved by 
using a Newton iteration. After the vector Xk converges, the ro- 
tation matrix Fk is obtained by (j52|) . 

This computational approach is desirable, since the implicit 
equations are solved numerically using operations in a linear vec- 
tor space. The three-dimensional rotation matrix Fk is computed 
by numerical iterations on , and its orthogonal structure is au- 
tomatically preserved by (j52|) . It has been shown that this compu- 
tational approach is so numerically efficient that the correspond- 
ing computational load is comparable to explicit integrators [15| . 
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= ^{jlAsljh' + EA)i\\q,4-hf, (46) 
= 0, (47) 
^Ml^Aqk^a-l + iMlAqk^a + ^Mf Z\gfc,,+i + i(M|i„ + M23^_i)Z\5p, - ^Mll^Aqk_^^^_^ - '^Ml_^Aqj,_^^^ 
-\Mll^Aqk-i,a+i - \{Mll^^^ + Ml\^_^)Asp^_, - -^(1 + 37V - 3a)Z^Sp, Z\gfc,«+i + -^Asp^Aq^^^ 

^hr>qk,a^k + ^(5 + SAT - 3a)Z\sp,Z\gfc,a-i = 0, (48) 
i(M2 + M)Z\gfc,iV+i + ^Mf Z\gfc,jv + ^Mf^Z^Sp, + ^MRk(Fk - I)pc - \{MI_^ + M)Aqk-i,N+l - ]- M^l^Aqk-i^N 
-^^fc-i,iV^^Pfc-i - \MRk_^{Fk_^ - I)pc + -^Asp^Aqk^N+i + ^Asp^Aqk,N + /iD,,^^^, = 0, (49) 
^(i^fc^d - JdF^ - JdFk-i + F^-iJdY + ^PcRl{Aqk,N+l - Aqk-i,N+l) - hMgpcRles = 0, (50) 



14. E. Hairer, C. Lubich, and G. Wanner, Geometric numerical 
integration, ser. Springer Series in Computational Mechanics 
31. Springer, 2000. 

15. T. Lee, M. Leok, and N. H. McClamroch, "Lie group varia- 
tional integrators for the full body problem in orbital me- 
chanics," Celestial Mechanics and Dynamical Astronomy, 
vol. 98, no. 2, pp. 121-144, June 2007. 

16. , "Lie group variational integrators for the full body 

problem," Computer Methods in Applied Mechanics and En- 
gineering, vol. 196, pp. 2907-2924, May 2007. 

17. T. Lee, "Computational geometric mechanics and control 
of rigid bodies," Ph.D. dissertation. University of Michigan, 
2008. 

18. T. Lee, M. Leok, and N. McClamroch, "Optimal attitude 
control of a rigid body using geometrically exact computa- 
tions on SO (3)," Journal of Dynamical and Control Systems, 
vol. 14, no. 4, pp. 465-487, 2008. 

19. B. Leimkuhler and S. Reich, Simulating Hamiltonian Dy- 
namics, ser. Cambridge Monographs on Applied and Com- 
putational Mathematics. Cambridge University Press, 2004, 
vol. 14. 

20. J. Marsden and M. West, "Discrete mechanics and varia- 
tional integrators," in Acta Numerica. Cambridge Univer- 
sity Press, 2001, vol. 10, pp. 317-514. 

21. A. Iserles, H. Munthe-Kaas, S. N0rsett, and A. Zanna, "Lie- 
group methods," in Acta Numerica. Cambridge University 
Press, 2000, vol. 9, pp. 215-365. 

22. L. Dieci, R. Russell, and E. van Vleck, "Unitary integra- 
tors and applications to continuous orthonormalization tech- 
niques," SIAM Journal on Numerical Analysis, vol. 31, no. 1, 
pp. 261-281, 1994. 

23. O. Junge, J. Marsden, and S. Ober-Blobaum, "Discrete me- 
chanics and optimal control," in IFAC Congress, Praha, 
2005. 



